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Abstract. Absorbing aerosols play an important, but uncer- 
tain, role in the global climate. Much of this uncertainty is 
due to a lack of adequate aerosol measurements. While great 
strides have been made in observational capability in the pre- 
vious years and decades, it has become increasingly apparent 
that this development must continue. Scanning polarimeters 
have been designed to help resolve this issue by making accu- 
rate, multi-spectral, multi-angle polarized observations. This 
work involves the use of the Research Scanning Polarimeter 
(RSP). The RSP was designed as the airborne prototype for 
the Aerosol Polarimetery Sensor (APS), which was due to 
be launched as part of the (ultimately failed) NASA Glory 
mission. Field observations with the RSP, however, have es- 
tablished that simultaneous retrievals of aerosol absorption 
and vertical distribution over bright land surfaces are quite 
uncertain. We test a merger of RSP and High Spectral Res- 
olution Lidar (HSRL) data with observations of boreal for- 
est fire smoke, collected during the Arctic Research of the 
Composition of the Troposphere from Aircraft and Satel- 
lites (ARCTAS). During ARCTAS, the RSP and HSRL in- 
struments were mounted on the same aircraft, and valida- 
tion data were provided by instruments on an aircraft fly- 
ing a coordinated flight pattern. We found that the lidar data 
did indeed improve aerosol retrievals using an optimal es- 
timation method, although not primarily because of the con- 
traints imposed on the aerosol vertical distribution. The more 
useful piece of information from the HSRL was the total 
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column aerosol optical depth, which was used to select the 
initial value (optimization starting point) of the aerosol num- 
ber concentration. When ground based sun photometer net- 
work climatologies of number concentration were used as an 
initial value, we found that roughly half of the retrievals had 
unrealistic sizes and imaginary indices, even though the re- 
trieved spectral optical depths agreed within uncertainties to 
independent observations. The convergence to an unrealis- 
tic local minimum by the optimal estimator is related to the 
relatively low sensitivity to particles smaller than 0.1 (pm) 
at large optical thicknesses. Thus, optimization algorithms 
used for operational aerosol retrievals of the fine mode size 
distribution, when the total optical depth is large, will require 
initial values generated from table look-ups that exclude un- 
realistic size/complex index mixtures. External constraints 
from lidar on initial values used in the optimal estimation 
methods will also be valuable in reducing the likelihood of 
obtaining spurious retrievals. 


1 Introduction 

Polarimetric remote sensing is a valuable tool for the study 

of atmospheric particles. In the 1970’s, Hansen and Hov- 

enier (1974) used observations of (linear) polarization in 

sunlight reflected from Venus to conclude that the planet 

is shrouded by sulfuric acid clouds (Hansen and Hovenier, 
1974). Today, the field has grown to include remote sensing 
of atmospheric aerosol particles surrounding our own planet. 
Aerosols are airborne particulate matter that are of interest 
because of the variety of ways they can affect the global 



Published by Copernicus Publications on behalf of the European Geosciences Union. 


7046 


K. Knobelspiesse et al.: Polarimeter and lidar combined retrievals during ARCTAS 


climate, and because this interaction is not well understood. 
Unlike the greenhouse gases associated with climate change, 
aerosols typically survive in the atmosphere for only days or 
weeks, and are thus highly regional. They have both nat- 
ural and anthropogenic sources and highly variable optical 
properties. For these reasons, the Intergovernmental Panel on 
Climate Change has identified aerosols as a major source of 
uncertainty in the radiative forcing of climate (IPCC, 2007). 
Much of this uncertainty is due to inadequacies in provid- 
ing appropriate and comprehensive information for climate 
models (Mishchenko et al., 2004, 2007b). 

The Aerosol Polarimetry Sensor (APS), which was part of 
the NASA Glory orbital mission, was intended to expand the 
quantity and quality of atmospheric aerosol descriptive pa- 
rameters. Unfortunately, the Glory spacecraft failed to reach 
orbit during its launch on 4 March 2011. APS was to use 
polarimetric observations of reflected sunlight at about 240 
viewing angles in each of nine channels at visible and near 
infra-red wavelengths (Mishchenko et al., 2007a). Aerosol 
property retrievals were to be determined by optimizing a ra- 
diative transfer simulation to match observed linear and po- 
larized reflectance. In order to prepare for the launch of APS, 
the Research Scanning Polarimeter (RSP) was constructed 
(Cairns, 2003). The RSP is similar to the APS, but is de- 
ployed on an airborne, rather than orbital, observation plat- 
form. The RSP has been used in a number of field campaigns 
in the last decade to test the ability to accurately measure 
a variety of aerosol types (Cairns et al., 1997; Chowdhary 
et al., 2001, 2002, 2005a,b; Cairns, 2003; Elias et al., 2004; 
Waquet et al., 2009). Despite the loss of Glory, RSP will 
continue to be deployed in the field. 

In June and July of 2008, the RSP participated in the 
summer phase of the Arctic Research of the Composition 
of the Troposphere from Aircraft and Satellites (ARCTAS) 
field campaign (Jacob et al., 2010). ARCTAS was the 
NASA contribution to the international Polar Study using 
Aircraft, Remote Sensing, Surface Measurements and Mod- 
els, of Climate, Chemistry, Aerosols, and Transport (PO- 
LARCAT) program of observations for the International Po- 
lar Year. While stationed in Northern Canada, the RSP 
(and other instruments associated with ARCTAS) observed 
smoke aerosols from boreal forest fires. Biomass Burn- 
ing (BB) aerosols tend to be both highly absorbing and un- 
evenly distributed vertically in the atmosphere. The effects 
of aerosol absorption and vertical distribution were found 
to have similar impacts on the polarized reflectances ob- 
served over bright desert surfaces which increases the un- 
certainty in the retrieval of both aspects of the aerosol dis- 
tribution (Waquet et al., 2009). Furthermore, retrievals of 
ocean optical properties may be affected by incorrect aerosol 
retrievals due to aerosol vertical distribution errors (Duforet 
et al., 2007). Orbital scanning polarimeter observations 
could therefore benefit from a data merger with active instru- 
ments such as the Cloud- Aerosol Lidar with Orthogonal Po- 
larization (CALIOP) onboard the Cloud-Aerosol Lidar and 


Infrared Pathfinder Satellite Observation (CALIPSO) satel- 
lite (Winker et al., 2003, 2004, 2007). CALIOP can deter- 
mine the vertical layering of aerosols and reduce the ambi- 
guity in BB aerosol retrievals due to aerosol vertical distri- 
bution uncertainty. During ARCTAS, the RSP flew onboard 
the NASA B-200 aircraft, which also carried the High Spec- 
tral Resolution Lidar (HSRL) (Hair et al., 2001, 2008, and 
Rogers et al., 2009). In addition, the B-200 flew coordi- 
nated missions with an aircraft containing instrumentation 
whose data can validate RSP results. ARCTAS was there- 
fore the ideal mission to evaluate combined polarimeter and 
lidar aerosol property retrievals for near source smoke, and 
that is the main goal of this paper. 

In the next section, we will provide an overview of the RSP 
retrieval approach. We will also briefly describe the HSRL 
and in situ data collection and processing, and provide refer- 
ences to more detailed discussions. Next, we will present the 
results of our tests. RSP retrievals were performed with and 
without HSRL data, and were compared to in situ data to test 
their success. Finally, implications of these tests for future 
RSP/HSRL measurements will be discussed. 


2 Method 

2.1 RSP instrument specifics 

The Research Scanning Polarimeter is an airborne prototype 
for the Aerosol Polarimetery Sensor. The main goals of both 
RSP and APS are to retrieve a complete suite of aerosol 
and cloud microphysical parameters from orbit (Mishchenko 
et al., 2004, 2007b). Both instruments have similar character- 
istics, but since the data we analyze in this paper were from 
RSP, that instrument alone will be discussed here. The RSP 
has nine optical channels with center wavelengths of 410, 
470, 555, 670, 865, 960, 1590, 1880 and 2250nm, and is a 
passive, along track, scanning (not imaging), device. Each 
RSP scan begins about 60° forward of nadir in the direction 
of aircraft motion, and samples at 0.8° intervals to about 60° 
aft of nadir. The instantaneous field of view (IFOV) of the 
RSP is fourteen milliradians. Lor an average (during ARC- 
TAS) aircraft height of 8680m above sea level and ground 
elevation of 370 m, this corresponds to ground pixel size of 
about 120 m. Lor observation of elevated smoke plumes such 
as ours, the spatial resolution is about 65 m. 

The RSP/APS observes linearly polarized and total re- 
flectance. This can be described by the first three terms of 
the Stokes polarization vector (Hansen and Travis, 1974), /, 
Q and U. I is the total radiance, while Q and U indicate 
the direction and magnitude of the linearly polarized radi- 
ance. The fourth term in the polarization vector, V, rep- 
resents circular polarization, which is very small for atmo- 
spheric aerosols (Kawata, 1978, for spherical particles), and 
is neglected in the analysis of RSP data. Prom a remote 
sensing standpoint, it is useful to express the polarization 
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components of the Stokes vector in terms of reflectance. This 
takes the exo-atmospheric radiance into account, and is cal- 
culated as follows. 


Ri 


l * rl 

F 0 cos 0 S 


r q = 


Q * i'o 
F 0 COSda 


0 ) 


Ru 


V * i 'o 

F 0 cos0 s 


where F 0 is the annual average solar exo-atmospheric irra- 
diance (W m -2 ), spectrally weighted for the particular RSP 
band, r 0 is the solar distance in AU (thus compensating for 
solar distance deviation from average throughout the year), 
and 9 S is the solar zenith angle. The RSP instrument has 
a high accuracy (0.2%) of Rq and Ry relative to Rj, as 
the same detectors are used to measure intensity and linear 
polarization, and because calibration of their relative gain is 
performed with each scan. 

For each scene, the RSP makes observations of R/, Rq 
and Ry in nine wavelengths for 150 view angles. This pro- 
vides upwards of four thousand measurements that are used 
to determine the aerosol optical properties that provide an 
optimal fit to the observations. In practice, far fewer mea- 
surements are used (for reasons described below), and there 
is also a strong correlation between observations at differ- 
ent view angles. However, this is sufficient information to 
retrieve the half dozen or so optically important parameters 
for each aerosol size mode. Furthermore, the broad spectral 
range of the RSP observations allows the surface to be con- 
strained almost independently of the atmospheric state. 


2.2 ARCTAS and data selection 


The ARCTAS field campaign involved three aircraft in two, 
three week segments in 2008. The RSP was deployed during 
the summer (June-July) phase of the campaign (ARCTAS-B) 
on a Beechcraft King Air B-200 aircraft based at the NASA 
Langley Research Center in Hampton, Virginia. One of the 
primary goals of ARCTAS-B was to observe and character- 
ize boreal forest fire smoke (Jacob et al., 2010), so the B- 
200 was based in Yellowknife, Northwest Territories, Canada 
(62° N, 1 14° W) for this purpose. The B-200 flew coordi- 
nated flights with another ARCTAS aircraft during satellite 
overpass times. This aircraft, a Lockheed P-3 Orion from 
the NASA Wallops Flight Facility in Wallops, Virginia, car- 
ried a variety of instrumentation, including the Hawaii Group 
for Environmental Aerosol Research (HiGEAR) aerosol in 
situ sampling instruments. Furthermore, the Ames Airborne 
Tracking Sunphotometer (AATS-14), a fourteen channel sun 
tracking sun photometer, was installed on the P-3. RSP data 
collected during these times is ideal for the validation of 
RSP aerosol retrievals of vertically inhomogeneous, highly 
absorbing aerosols. We chose to study a small data subset 
from the smoke plume of a recent boreal forest fire. Surveys 


of other boreal forest fire aerosols can be found elsewhere 
(e.g. Eck et al., 2003, 2009; Koppmann et al., 2005; Reid 
et al., 2005a,b). 

A scene on 30 June 2008, was selected for analysis. On 
that day, the B-200 overflew a smoke plume downwind from 
its source in conditions that were otherwise nearly free of 
aerosols. The fire was in northern Sasketchewan, Canada 
(58.41° N, 106.81° W) and had been burning since 29 June. 
Natural Resources Canada, in its forest inventory of 2001 
(CanFI2001), characterizes this region as a mixed soft and 
hardwood boreal forest. RSP observations were made about 
130 km downwind from the source (see Fig. 1). Back trajec- 
tory analysis (Fuelberg et al., 2010) indicates that the smoke 
was 2h and 17 min old at this point, although it should be 
noted that this analysis did not account for local dynamics 
associated with the fire. We consider this scene to be ideal 
for a number of reasons, as described below. 

1 . The HSRL was operational at the time of RSP observa- 
tion, indicating an optically thick plume detached from 
the surface. 

2. The B-200 aircraft (containing the RSP and HSRL) was 
in a coordinated flight path with the P-3 aircraft (con- 
taining AATS and HiGEAR), which was under-flying 
the B-200 at a mean altitude of 627ma.s.l. (above sea 
level). 

3. A MODIS-Aqua image was collected about forty min- 
utes prior to the scene. The “truecolor” version of this 
scene provides the spatial context as observed in Fig. 1. 

4. Aircraft orientation was close to the solar principal 
plane (the plane containing both the solar illumination 
and observation vectors), providing a large range of 
scattering angles for use during optimization of RSP 
data. 

5. The atmospheric state outside the plume was excep- 
tionally free of aerosols. An RSP aerosol retrieval 
performed prior to contact with the plume yielded an 
aerosol optical thickness of 0.07 at 555 nm, while AATS 
observations from the same location found an optical 
thickness of 0.04 at 520 nm. 

6. Aircraft geometry, as provided by an onboard inertial 
monitoring unit, had relatively low uncertainties. 

7. Back trajectory analysis (Fuelberg et al., 2010), com- 
bined with the Canadian Wildland Fire Information Sys- 
tem database (Canada, 2009) shows that the airmass in 
the days prior to the observation did not come into con- 
tact with any significant sources of biomass burning (or 
other) aerosols. 
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Fig. 1. MODIS Aqua satellite fire imagery from 19:50 UTC. The RSP observations, indicated by the ordered white polygons at the right, 
were collected at 20:28 UTC on the B-200 aircraft. The flight track of the P-3 aircraft is shown in red. The portion of the P-3 flight under the 
indicated RSP observations was coordinated to occur at the same time. Google Earth was used to visualize these data. 


Aerosol retrieval optimization was performed individually 
for about twenty segments in each scene. To reduce noise, 
each segment is an average of five RSP scans. We used this 
approach (rather than performing a single optimization on an 
average of a larger number of scans) to test the consistency 
of the optimization approach, and to allow for the compensa- 
tion, or quantification, of the effects of scene heterogeneity 
and changing aircraft geometry. In situ and sun photometer 
data on the P-3 indicate that aerosol optical thickness varies 
across the smoke plume, while intensive properties such as 
single scattering albedo remain mostly constant (see figures 
in Sect. 3). We would therefore expect retrieved aerosol 
number concentration, n, (and thus optical thickness, r), to 
vary across the plume, whereas intensive properties (size and 
composition) should vary less. For the time, location, solar 
geometry and other details associated with each segment, see 
Table Al. 

2.3 Optimization approach 

We retrieved aerosol optical properties by matching the out- 
put from a radiative transfer model to polarized and unpo- 
larized reflectances observed by the RSP. The Levenberg- 
Marquardt approach, which is a standard method for the 
solution of non-linear least squares minimization problems, 
was used to find that match. With the exception of some 
methodological and software differences (outlined below), 
we used the approach described in detail in Waquet et al. 
(2009). Generally speaking, the radiative transfer model and 


software is identical, using Lorenz-Mie computations of sin- 
gle scattering properties and the Doubling and Adding ap- 
proach to create multiple scattering in a layered atmosphere 
(Mishchenko and Travis, 2008; Hansen and Travis, 1974). 
Optimization is performed differently than in Waquet et al. 
(2009), with a slightly different cost function and measure- 
ment vector and using publicly available software to perform 
the Levenberg-Marquardt optimization. This software, im- 
plemented in the Interactive Data Language (IDL) as MP- 
FIT (Markwardt, 2009), is a translation of the original FOR- 
TRAN language MINPACK-1 (More, 1977) solver, and is 
well documented. For this reason, we chose to use MPFIT, 
as it is faster and more flexible than the optimization soft- 
ware created for Waquet et al. (2009) (although both tech- 
niques are appropriate). Furthermore, numerical and compu- 
tational issues associated with techniques such as this have 
been addressed in the creation of MINPACK-1 and MPFIT. 
However, there are often subtle differences between different 
implementations of the Levenberg-Marquardt approach, we 
describe our methodology in more detail below. 

2.3.1 Optimal estimation 

Optimal estimation is performed by the minimization of a, 
generally nonlinear, cost function. The cost function incor- 
porates information about the uncertainty in the measure- 
ments so that the estimate that minimizes the cost function 
is optimal in the sense that it provides the best possible fit 
to the measurements given their uncertainties. In our case. 
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this cost function, <t>, is the sum of squares of differences be- 
tween the observational data and model calculation vectors 
weighted by their total uncertainty. 


<*>(*) 


\ II WII 2 


1 

2 


Y - G(x) 2 
Cj 


( 2 ) 


where Y is the measurement vector and G(A)is the modeled 
vector for aerosol parameters x. We use ||.|| to indicate the 
use of the Euclidean norm. Note that, unlike Waquet et al. 
(2009), this cost function depends only on the differences be- 
tween observational data and model results, and does not use 
and additional side constraint on the distance from the ini- 
tial model parameters. In this sense, it is free to explore the 
parameter space available during optimization with a mini- 
mized dependence on initial conditions. Cj is the measure- 
ment error covariance matrix, which can be broken down into 
component errors as 


c T — C £ + Ccal + Cpol + C ag + Cy + Cp. (3) 

Following the example of Waquet et al. (2009), the in- 
strumental noise, C £ is shot noise limited, C ca i is the ab- 
solute radiometric calibration uncertainty of 3 %, and C po ] 
accounts for the uncertainty in the relative polarimetric accu- 
racy, which increases with the degree of linear polarization. 

C e>/ = 10 -7 cosP s Rf(i) 


Ccai,,- = (0.03 Rq( i)) 2 (4) 

Cpol,/ = (0.001 [R,(i) + |/? G (i)|]) 2 

where the index, i, indicates that all the uncertainties are 
spectral band and view angle dependent and it is assumed 
that the measurement error covariance matrix is diagonal. 

Waquet et al. (2009) uses a fourth term in Eq. (3) which 
accounts for errors modeling the surface reflectance. They 
label this term as Cf- The contribution of Cf to the total er- 
ror is very small, and since the smoke plume we observe is 
quite optically thick, for this case we can (and do) neglect 
Cf- For the data we analyze here, we found two additional 
sources of uncertainty that must be included in the error bud- 
get for the observations. The variability encountered when 
averaging several scans of an inhomogeneous plume requires 
a new term,C ag , which represents the sampling uncertainty 
for that observation and is the standard deviation for the av- 
eraged set of scans for each view zenith angle. We have also 
added two terms to account for uncertainties in the aircraft 
geometry. C y is the uncertainty in the observed polarized re- 
flectance caused by uncertainties in the aircraft heading angle 
(yaw), which expresses itself as an incorrect view azimuth 
angle. We estimated the magnitude of this uncertainty by 
shifting the observed data for a heading error of 0.5°, and 
taking the difference between the observed and shifted ob- 
servation. The 0.5° shift was empirically determined based 
on data variability. For Rq, we add to C y the error from an 


incorrect rotation into the scattering plane. C p is the uncer- 
tainty in the observed polarized reflectance caused by aircraft 
pitch angle, which expresses itself as an incorrect view zenith 
angle. Uncertainty is computed in a manner similar to C y , 
and we again used a value of 0.5°. Both estimates of geom- 
etry uncertainty were determined empirically based on data 
variance. Typically, the largest contributer to measurement 
uncertainty is aircraft pitch, closely followed by sampling 
uncertainty due to scene heterogeneity. Absolute calibration 
uncertainties are usually a much smaller contributer to the 
observation errors than the uncertainties in aircraft geometry 
and scene variability. 

The above covariance matrix definition applies to com- 
ponents of the measurement vector associated with the po- 
larized reflectance. The covariance for the total reflectance 
is slightly different. Polarimetric accuracy does not affect 
the uncertainty in the total reflectance, instrumental noise re- 
mains the same, and the radiometric calibration uncertainty 
is now 3 % of the total reflectance. 

C cal (0 = [0.03 R/(t)] 2 (5) 


The assumption of uncorrelated measurement errors is 
commonly used (Lebsock et al., 2007; Hasekamp, 2010) and 
is valid for the instrument noise, and to a lesser extent for 
polarimetric accuracy. For radiometric accuracy and aircraft 
attitude, there will be correlation between measurement un- 
certainties at different view angles. For the data analyzed 
here, however, scene heterogeneity and random variations in 
aircraft attitude dominate the error budget and justify the as- 
sumption of uncorrelated measurements errors. Future work 
will incorporate error correlation for aerosol retrievals over 
the ocean, where a lack of surface heterogeneity does not 
overshadow the effect of correlations between measurement 
errors for different view angles. 

Starting from selected initial values for the aerosol param- 
eters, x 0 (see Sect. 2.3.3), the Levenberg-Marquardt algo- 
rithm provides an iterative search for the parameter vector, 
p, that minimizes the cost function <J> (x 0 + p) and thus gives 
the optimal solution x = (jc 0 + p)- The iteration proceeds by 
constraining the minimization of the cost function at any step 
to lie within a trust region. A*, around the current iteration 
step. A linearization about the current state vector, Xk, is 
used to determine the size of the next step. At any given 
step, we are therefore attempting to solve 

min {||F(x*) + J* p k \\ : ||D* p k \\ < A a .} (6) 


where k is an index indicating how many steps have been 
made, and J k is the Jacobian matrix, which expresses the 
forward model sensitivity to parameter change. The notation 
: means “such that”. 


J* 


3F(x) 


dx 


x=x k 


(7) 


D k is a diagonal, but not an identity, matrix. This was in- 
troduced by Marquardt to allow for a large step in a direction 
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with low curvature and a small step in a direction with high 
curvature. Heuristically, this serves to reduce the effects of a 
narrow valley with a relatively flat floor and is calculated as 
part of the iteration with the formula 


Da = diag (d[ k) , ..., 4*°) 


,(A) 

a) = max 


1 

1 

3F(x) 

dxi 


L 


x=x k 


( 8 ) 


The Levenberg-Marquardt method is based on the theorem 
that if p* is a solution to Eq. (6), then p*4 — p( T) for some 
T > 0 where 


pm = - (J l Ja + T D[ D,)- 1 J l F(x k ) (9) 


The constrained minimization is then implemented by the 
following algorithm, from More (1977) 


1. Given A k > 0, find T k > 0 such that if 


(J 2 Ja- + Xt-Dj Dyfc) p k = -J l F (x k ) (10) 


This iteration continues until either the step size is less 
than the expected uncertainty in state (aerosol parameter) 
space, or the reduction in the cost function is small. Once 
the iteration is complete, the retrieval error covariance ma- 
trix, C A , can be obtained from the Jacobians computed in the 
final step using the equation 

C, = (J r C r J)” 1 . (12) 

Uncertainty associated with retrieved parameters is the 
square root of the diagonal elements of C x . The error co- 
variance matrix can also be used to compute retrieval uncer- 
tainty of a parameter that is not directly retrieved during opti- 
mization (such as the aerosol optical thickness, see Table 2), 
provided the dependent variables of the required parameter 
are part of the retrieval vector x. This is done by computing 
the sensitivity of the indirectly retrieved parameter in the for- 
ward model to changes in each of the retrieved parameters. 
Following the method of Hasekamp and Landgraf (2007), 
the uncertainty in the indirectly retrieved parameter (denoted 
with the subscript A) is 


then either T k =0 and IID^/^H < A k , or Y,t>0 and 
(1 — cr)A k < lIDjp^H < (1 +er)A k . o is the relative er- 
ror within which the iterative estimate of T k is required 
to meet HD^H = A k . 

2. If ||F(oca + p k )\\ < ||F(*a)|| then set x k+ \ =x k + p k . 
Otherwise, shrink the trust region ( A k+ 1 = j A k ) and re- 
turn to step 1, without modifying x and J. 

3. Choose A^+j. The ratio of the actual reduction of the 
cost function to the predicted reduction i s the criterion 
for adjusting A, viz., 

= ||F(a:a)1| 2 - \\F(x k + p k ) || 2 

Pk ||F(x *)|| 2 - ||F(x*) + JaP, II 2 ' ' ' 

The value of p k is used to determine A k+ \ as follows 

- if p k < 1/4, then set Aj.+i < A k using the method 
in Fletcher (1971). 

- if \ < Pk < \ and T k = 0, set A k +\ =2\\V k p k || 

- if \ < Pk < l and T k > 0, set A k+l = A k 

- if Pk > 3/4, then set A^+i = 2 1 1 p k ||. 

To summarize, if the predicted and actual reductions in 
cost functions are similar then the trust region is in- 
creased, whereas if they are widely disparate the trust 
region is decreased based on a prescription by Fletcher 
(1971). 

4. Update D^+i using Eq. (8). 


N N 


E E c '.u/ 


dG a x) 

dxj 


(13) 


where i and j are subscripts identifying elements of the re- 
trieved parameter vector (which has a length of N). G a (x ) is 
the forward model value for the indirectly retrieved parame- 
ter, A, evaluated at the optimal estimate, x (note that this is a 
scalar value not in the same space as G(jc)). We compute the 
forward model sensitivity for an indirectly retrieved param- 
eter ( ^%-— ) numerically, with a post-processing step after 
optimization is finished. 

We chose this form of the Fevenberg-Marquardt method 
for several reasons, which we briefly describe here. For 
more details about this method and its convergence prop- 
erties, see More (1977), or the documentation of the MIN- 
PACK FORTRAN subroutines on which our technique is 
based (www.netlib.org/minpack/ and More et al., 1980). The 
Fevenberg-Marquardt technique is robust, and is efficient in 
its use of forward model calculations. Our vertically inho- 
mogeneous vector radiative transfer model is computation- 
ally intensive, so it is not feasible to make multiple forward 
calculations to determine an optimal parameter T k . The al- 
gorithm given above uses an iterative estimate of T k , but the 
iteration is one that finds a step within a trust region that 
does not require additional forward model or Jacobian calcu- 
lations (Hebden, 1973). For a trust region criterion of 10% 
(A = 0.1), T is typically found within two iterations (More, 
1977). The trust region grows whenever the predicted and 
actual cost function reductions are similar, so the transition 
between the steepest descent and Newton-Gauss components 
of the Fevenberg-Marquardt technique is heuristically rea- 
sonable, as whenever the next step is within the trust region 
an efficient Newton-Gauss step is used. 
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2.3.2 The measurement vector 

The measurement vector, F, contains observations of the Q 
component of polarized reflectance (Rq) at six wavelength 
bands, and total reflectance ( Rj ) for one band. In each band, 
observations were taken from a scan between 20° forward in 
the aircraft and 40° aft. Since the RSP has an angular res- 
olution of 14milliradians, this corresponds to 75 individual 
measurements in each scan, for a total of 450 observations of 
Rq and 75 observation of Rj. The polarization frame of ref- 
erence was rotated into the scattering plane (the plane con- 
taining both the solar and observation vectors) so that Rq 
measurements alone are sufficient to represent the linear po- 
larization (Rjj is negligible for single scattering in the scat- 
tering plane). For the purposes of optimization using data 
near solar backscatter, the choice of Rq is preferable to the 

polarized reflectance (R p = J R 2 q + Ry) that was used in Wa- 
quet et al. (2009). Rq is a signed measure and makes more 
effective use of observed neutral points than R p . 

The six wavelength bands used for Rq were centered at 
410, 470, 555, 670, 865 and 1590 nm. Three RSP bands were 
not used. Bands centered at 960 and 1880 nm are designed 
to allow the estimation of column water vapor and the iden- 
tification and characterization of cirrus clouds, respectively. 
The 2250 nm band is used to characterize the total and polar- 
ized surface reflectance for the radiative transfer model (see 
Sect. 2.3.3). Our particular scene contains optically thick 
smoke aerosols, so we also used one short wavelength Rj 
band (410 nm) in the optimization. Observations at this band 
have almost no contribution from surface reflectance, since 
surface reflectance is minimal in the blue unless it is covered 
by snow or ice (Hsu et al., 2006; Guanter et al., 2008) and 
the large aerosol load obscures what is reflected. Neverthe- 
less, we used the methodology of Kaufman et al. (1997) to 
characterize R / at 410 nm based on observations at 2250 nm, 
but should be insensitive to the uncertainty in this approach. 

The measurement vectors were averages of five consec- 
utive scans (a “segment” as described above). Each scan 
is composed of data that have been reorganized so that all 
views point to a specific location, at a specific altitude, in 
the atmosphere. This is essential to reduce heterogeneity in 
the aerosol properties and loading within each segment. We 
performed this data reorganization so that each segment con- 
tains all of the observations, at different zenith angles, of a 
point at an altitude of 3800 m. 

2.3.3 The radiative transfer model 

The radiative transfer model, G(a:), is actually two nested 
models. The inner model computes single scattering prop- 
erties of bi-modal log-normal size distributions of aerosols 
using a Mie code for spheres (Hansen and Travis, 1974). 
The elements of the retrieval vector required by the sin- 
gle scattering model are the complex refractive index, m. 


effective radius, r e , and effective variance, v e , for each size 
mode. The smaller size mode, which we denote hereafter 
with a f subscript, represents “fine” aerosols that are typ- 
ically the product of chemical processes, while the larger, 
“coarse” mode, denoted hereafter with a c subscript, repre- 
sents larger aerosols that are the result of a mechanical pro- 
cess. We expect the fine mode to dominate in our scene, 
since smoke aerosols are chemically generated. We also ex- 
pect that Mie models of spheres are adequate expressions of 
our aerosols. Indeed, the HSRL depolarization observations 
are very low, indicating spherical particle dominance (see 
Sect. 2.4.1). Furthermore, it is believed that smoke aerosols 
that have aged for an hour or more generally collapse into 
compact shapes that are well described, optically, by spheres 
(Martins et al., 1998; Reid and Hobbs, 1998; Abel et al., 
2003; Liu et al., 2008). The doubling and adding (DA) tech- 
nique is then used to compute the upwelling, downwelling 
and reflected polarized radiance fields generated by multi- 
ple scattering in an inhomogeneous atmosphere (Hansen and 
Travis, 1974; De Haan et al., 1987). The additional retrieval 
vector elements required for DA calculation are the vertical 
distribution of aerosol number concentration for each mode 
(/if) and ( n c ). The result is assembled into a vector represent- 
ing the geometry and wavelengths required to fill G(jc), and 
used in the optimization. 

The radiative transfer model also requires that the ground 
reflectance is specified. For polarized reflectance, we can 
utilize established assumptions about surface polarized re- 
flectance. Polarized reflectance of natural surfaces (such as 
the epicuticular wax coating the surfaces of leaves or the min- 
eral facets of exposed rock) will be dominated by interactions 
at the surface, as bulk scattering within the material tends to 
be weakly polarized. Fresnel reflectance off surface facets 
is thus an appropriate basis for modeling the polarized sur- 
face reflectance (Vanderbilt et al., 1985; Grant, 1987; Van- 
derbilt and Grant, 1991; Breon et al., 1995; Nadal and Breon, 
1999; Cairns, 2003; Elias et al., 2004; Waquet et al., 2009). 
Since Fresnel polarized reflectance depends mainly on ge- 
ometry and the real component of the refractive index, it ex- 
hibits little spectral variance. This is because most surface 
materials have a minimal real refractive index spectral vari- 
ance (the imaginary component of the refractive index does 
vary spectrally and gives total surface reflectances color). We 
can therefore use the longest wavelength channel (2250 nm), 
where aerosol effects are smallest, to characterize the sur- 
face reflectance at all wavelengths. This is used as the lower 
boundary condition in the DA radiative transfer model. In 
practice, the Fresnel polarized reflectance model is scaled to 
an appropriate value to match observed reflectances for each 
scene in order to account for surface roughness or variability 
in the real refractive index. We use a single scaling coef- 
ficient for all geometries, which differs from Waquet et al. 
(2009), who used a value that varied with geometry. This 
simplification is feasible for our analysis because of the rela- 
tively weak contribution of the polarized surface reflectance 
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due to optically thick aerosols, and allows the use a data 
closer to the solar backscattering angle. 

For unpolarized ground reflectance, observations at 
2250 nm were fit to the “RossThick-LiSparse” kernel 
based, bidirectional reflectance distribution function (BRDF) 
model. Fitting to the BRDF model mimicked the NASA 
Moderate-Resolution Imaging Spectroradiometer (MODIS) 
surface reflectance retrieval methodology (Lucht et al., 2000) 
and was performed previously with the RSP instrument in 
central Oklahoma in the United States (Knobelspiesse et al., 
2008). The results were multiplied by 0.25 and used as 
the reflectance at 410 nm. This is similar to the approxi- 
mation of Kaufman et al. (1997) for MODIS retrievals of 
aerosol properties over land. While it has been found that 
the 0.25 scaling is not appropriate for all surfaces and view- 
ing angles (Gatebe et al., 2001 and Remer et al., 2001), the 
variation is small enough to be rendered irrelevant by the op- 
tically thick aerosol layer obscuring the surface. For exam- 
ple, the retrieved aerosol optical thickness at 410 nm is 1.03 
for the sample optimization in Fig. 2. This corresponds to 
a mere 36% transmittance at nadir, not including the ef- 
fect of molecular scattering. Reflectances at 2250nm were 
typically low, averaging about 5 %. According to Kaufman 
et al. (1997), the 410 nm reflectance would be one quarter of 
this, 1.25 %. This surface reflectance transmitted through the 
aerosol layer at nadir would only be about 0.5 in reflectance 
units, while the total unpolarized reflectance is typically be- 
tween 15 % and 20,%. Errors associated with the approxi- 
mation of 410 nm reflectance would thus be suppressed in a 
similar manner. For this reason, one short wavelength band 
was included in the optimization vector G(;c). Longer wave- 
length bands were not used because of the risk of sensitivity 
to inaccurate surface reflectance characterization, both due to 
higher surface reflectance values in the red and near-infrared, 
and lower surface-obscuring aerosol transmission. 

Table 1 lists the model parameters, their units, and the a 
priori optimization values (x 0 ). A priori values were selected 
using the boreal forest fire smoke properties in Dubovik 
et al. (2002), which were derived from systematic obser- 
vations of sunphotometers in the Aerosol Robotic Network 
(AERONET). Parameters listed in italics in Table 1 are re- 
trieved during optimization, while the others are fixed. 

It has been suggested that the imaginary refractive index 
of some types of biomass burning smoke have a spectral de- 
pendence, specifically an increase in absorption (imaginary 
refractive index) in the blue and ultra-violet (Andreae and 
Gelencser, 2006; Lewis et al., 2008; Russell et al., 2010; 
Bergstrom et al., 2010). Unfortunately, direct observations 
of smoke complex refractive indices are very limited. Ab- 
sorption in smoke aerosols is typically associated with Black 
Carbon (BC) whose imaginary refractive index is spectrally 
flat. Organic Carbon (OC), which covers a wide range of 
species with varying chemical concentrations, is the source 
of absorption at short wavelengths (Kirchstetter et al., 2004). 
The relative contribution of BC and OC varies considerably 


Table 1. Retrieved aerosol parameters and the initial values used 
in optimization. Parameters listed in italics are optimized during 
the retrieval, while all others are kept constant. a Refractive In- 
dex values have no spectral dependence. b For optimizations that 
utilize data from the HSRL, the a priori number concentration is 
determined by the HSRL observations. c For optimizations that uti- 
lize data from the HSRL. the aerosol layer height is fixed by HSRL 
observations and not changed during optimization. 


Parameter 

Symbol 

Unit 

a priori 

Fine mode parameters 

Real refractive index 3 

1Z (wf) 

nidi 

1.52 

Imaginary refractive index 3 

1 (Wf) 

nln 

0.0094 

Effective radius 

r e,f 

(am 

0.15 

Effective variance 

v e.f 

nidi 

0.20 

Number concentration 13 

n f 

—2 

mum 

5.5 

Coarse mode parameters 

Real refractive index 3 

1Z (w c ) 

nidi 

1.52 

Imaginary refractive index 3 

1 (m c ) 

nidi 

0.0094 

Effective radius 

r e,c 

(am 

3.21 

Effective variance 

^e,c 

n/a 

0.23 

Number concentration 13 

n c 

_2 

mum 

0.0001 

Other parameters 

Aerosol height 3 

h 

m 

5000 


depending on the fire type and smoke age. BC is usually 
more associated with hot, flaming fires, while OC is greater 
in cooler, smoldering fires. Since our smoke plume was 
created by young fires with many hotspots, we expect that 
BC will dominate. Initially we did include a retrieval pa- 
rameter that allowed the refractive index to increase linearly 
(as wavelength decreases) for wavelengths less than 532 nm. 
However, most of our retrievals converged to a solution with 
no imaginary refractive index spectral dependence. Seg- 
ments that did converge with a spectral dependence typically 
had a large residual error. We therefore conclude that ei- 
ther an imaginary refractive index spectral dependence does 
not exist for this smoke, or we are not sensitive to it. We 
therefore used a spectrally flat imaginary refractive index in 
our retrievals. 

To facilitate comparison of our results with data collected 
by other instruments, we computed a variety of other aerosol 
properties using the retrieved aerosol properties. These de- 
rived parameters are listed in Table 2, along with information 
about the method used to derive them if they are not byprod- 
ucts of the forward model. Retrieved parameter uncertainty 
was determined from Eq. (12) as described above, and de- 
rived parameter uncertainty was calculated numerically as in 
Eq. (13). 
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Table 2. Derived aerosol parameters. A is wavelength. 


Parameter 

Symbol 

Unit 

Calculation 

Scattering cross-section 

Os (A) 

2 

(jrrr 

model byproduct 

Absorption cross-section 

o- a (A) 

2 

prrr 

model byproduct 

Extinction cross-section 

o'e (A) 

2 

(jm 

ffe (A) = ff a (A) + os (A) 

Asymmetry parameter 

g (A) 

n/a 

model byproduct 

Aerosol optical thickness 

r (A) 

n/a 

r (A) = n ff e (A) 

Angstrom exponent 

a 

n/a 

Slope of In (A) vs. In (er e (A)) 

Single scattering albedo 

m (A) 

n/a 

m (A) = a s (A )/cr e (A) 

Phase function at backscatter 

P n (A, 0 = 180°) 

n/a 

model byproduct 

Backscatter extinction ratio 

5(A) 

sr 

4 n/[nr (A) P n (A, 0 = 180°)] 


2.3.4 Atmospheric layer heights 

As discussed in Waquet et al. (2009), absorbing aerosol prop- 
erties are difficult to retrieve without information about their 
vertical distribution. This paper is therefore intended to test 
how to best model absorbing aerosol vertical distribution. To 
do so, we performed optimizations with and without data 
supplied by the HSRL, and compare the results. 

Retrievals performed without HSRL data modeled 
aerosols in a single, uniform layer attached to the ground. 
The top of the layer was allowed to vary as an optimiza- 
tion parameter. As listed in Table 1, the initial altitude was 
5000 m, while initial values for aerosol number concentra- 
tions were selected to replicate the mean optical depth for 
biomass burning aerosols as given by Dubovik et al. (2002). 

Unlike traditional backscatter lidars, the HSRL has 
the ability to discriminate between aerosol and molecular 
backscatter in the 532nm channel. This provides an inde- 
pendent estimation of aerosol extinction and backscatter co- 
efficients, and means that the observed aerosol optical depth 
is directly estimated with no microphysical assumptions re- 
quired. A more detailed description of this instrument can 
be found in Sect. 2.4.1. To find aerosol layer top and bot- 
tom heights, we took the derivative of the aerosol volume 
backscatter coefficient, b, after it had been weighted by the 
two pass atmospheric transmission (using the extinction co- 
efficient, k e ) from the specified altitude, h, to the aircraft alti- 
tude, h a (Eq. 14). Arbitrary thresholds were set for the result, 
y(h), where absolute values above the threshold indicate an 
aerosol layer top or bottom. 

y(h) = — (b{h) (14) 

d h V / 

Aerosol layer heights were fixed within the optimization, 
and aerosols were distributed evenly with respect to pressure 
throughout each layer. The results of this layer identification 
can be seen in the first panel of Fig. 3. Red dashed lines in- 
dicate the tops and bottoms of the pair of layers identified 
using Eq. 14. An additional piece of information from the 


HSRL is the total column aerosol optical depth at 532 nm. 
This was used to select a priori n c and n f, which were dis- 
tributed according to the aerosol optical depth in each layer. 
The relationship between n and aerosol optical depth was 
determined by the extinction cross-section calculated by the 
forward model for initial values of refractive index and size. 
Ratios of aerosol optical depth for both modes between dif- 
ferent layers were maintained throughout the optimization. 

2.4 Validation data 

Results of the optimization with and without HSRL data 
were compared to a variety of contemporaneous measure- 
ments by other instruments. The first comparison we made 
was between optimization results and HSRL observations of 
aerosol optical thickness, r, and the backscatter to extinc- 
tion ratio, S, at 532 nm. The former comparison is somewhat 
compromised by the fact that it was used to set the initial 
aerosol number concentration. However, since the initial val- 
ues for size and complex refractive index are different than 
the retrieved values, the optical depth of a specified aerosol 
number concentration must also change. All of these values 
are unconstrained, so agreement between RSP and HSRL 
optical depth is by no means guaranteed by the choice of 
the initial value. In addition to comparisons with the HSRL, 
the B-200 flew in coordination with the NASA P-3 aircraft, 
which carried two instruments of interest to us. The AATS- 
14 instrument provided the column aerosol optical thickness 
(above the aircraft) in fourteen narrow wavelength channels. 
The aerosol Angstrom exponent, a, was also calculated us- 
ing these data by performing a linear regression in the log-log 
domain. The HiGEAR suite of in situ sampling instruments 
provided the aerosol absorption (k a ) and scattering (k s ) co- 
efficients, from which the single scattering albedo, m is de- 
rived. HiGEAR instrument suite also measured aerosol size 
information, which we present as the effective radius, r e and 
effective variance, v e . Details about the collection of these 
data and the motivation for their use to validate RSP aerosol 
property retrievals are described below. 
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2.4.1 HSRL 

The HSRL, as described above and in more detail in Hair 
et al. (2001, 2008) and Rogers et al. (2009) is a dual band 
(532 and 1064 nm) lidar that was deployed along with the 
RSP on the B-200 aircraft. The HSRL uses an iodine 
cell with a narrow spectral absorption feature centered on 
the laser line at 532 nm to absorb radiation scattered by 
aerosols (close to the laser line) and transmit light scattered 
by molecules (Doppler shifted away from the laser line). 
Another channel is used to measure the light scattered by 
aerosols and molecules at 532 nm, similar to the usual elastic 
backscatter lidar measurement. This capability means that 
the HSRL can distinguish molecular from aerosol backscat- 
ter without the microphysical assumptions that are required 
for regular lidars. At 1064 nm, the typical lidar approach is 
used, although the 532 nm measurements can be used to con- 
strain the derivation of backscatter profiles (Hair et al., 2008). 
The HSRL also observes the depolarization ratio, d, which 
is the ratio of the perpendicular (cross-polarized) to paral- 
lel polarized aerosol backscatter coefficient. This parameter 
provides a measure of the sphericity of the aerosols since is 
a unitless ratio. In the case of single scattering it is zero for 
spheres. During this scene, the HSRL observed particulate 
depolarization of about 0.05, with is small enough that the 
aerosols can be treated as spheres (as they are in our radia- 
tive transfer model). The spherical assumption is reasonable 
for slightly aged smoke particles, which start as filaments and 
collapse into more compact shapes in the first few hours after 
creation (Abel et al., 2003; Reid and Hobbs, 1998; Liu et al., 
2008). 

2.4.2 AATS-14 

The Ames Airborne Tracking Sunphotometer (AATS-14) is 
a fourteen channel (354 nm to 2139 nm) autonomous sun 
tracking sun photometer. It is mounted externally on the 
top of an aircraft, and provides above aircraft aerosol opti- 
cal thickness by continuously tracking the direct solar beam. 
This is successful provided the solar disk is not obscured by 
clouds or interference from the aircraft body and that the air- 
craft angular motion does not exceed 6° per second. Cali- 
bration is extremely stable, typically maintained within 1 % 
over the course of about a year. The AATS-14 has been de- 
ployed on a variety of aircraft during field campaigns since 
1996 (for example, Russell et al., 1999; Schmid et al., 2003, 
2006; Redemann et al., 2005; Livingston et al., 2009). Dur- 
ing ARCTAS, the AATS-14 was mounted on the P-3 aircraft, 
which flew under the B-200 at a low altitude (Shinozuka 
et al., 2011). AATS-14 should provide an excellent mea- 
sure of the spectral dependence of aerosol optical depth. As 
we will see below (with the in situ instrument data), there 
are smoke related aerosols at the P-3 altitude, even if the to- 
tal quantity of aerosols between the P-3 and the ground is 
much less than the quantity above. Thus, the column aerosol 


optical thickness that we observe with the RSP is at least as 
large as the values measured by the AATS-14, and should be 
slightly higher. 

For validation purposes, we want to compare aerosol opti- 
cal thickness at the same wavelength. The RSP optimization 
software can produce an optical thickness at any wavelength 
(with an extra radiative transfer model run with the final pa- 
rameters at the specified wavelength), while the HSRL only 
provides an independent measure of column aerosol optical 
thickness at 532 nm. This is not an AATS channel (the clos- 
est is 520 nm), so AATS aerosol optical thickness at 532 nm 
was determined by fitting a second order polynomial to the 
spectra in the log-log domain. Uncertainty values were inter- 
polated in the same manner. 

2.4.3 HiGEAR 


HiGEAR is a suite of instrumentation deployed by the Uni- 
versity of Hawaii on the P-3. HiGEAR provided valida- 
tion data in the form of absorbing and scattering coefficients 
(from which the Angstrom exponent and single scattering 
albedo can be derived) and the aerosol particle size distri- 
bution. Some examples of results from previous field cam- 
paigns can be found in Clarke et al. (2007) and McNaughton 
et al. (2009). Shinozuka et al. (2011) integrated the in situ 
measurements over the vertical profiles flown during ARC- 
TAS, and demonstrated that the resulting layer aerosol opti- 
cal depth agrees with AATS-14 measurements within 3 % or 
0 . 02 . 

The aerosol size distribution was determined by two 
instruments. The Long Differential Mobility Analyzer 
(LDMA) counts particles in the 0.01 to 0.5 qm diameter 
range (Clarke et al., 1998). The LDMA maintains the 
aerosols at ambient temperature and pressure, at dry (less 
than 30%) humidities. Particles in the 0.5 to 10 qm range 
were measured by a TSI model 3321 aerodynamic particle 
sizer (APS) (McNaughton et al., 2007). Like the LDMA, the 
APS operates at ambient temperature and pressure at dry hu- 
midities. Effective radius and variance were calculated by 
fitting a log-normal distribution to the data and deriving the 
effective radius and variance values from the fit parameters 
(see Eq. 2.60 and the next page of discussion in Hansen and 
Travis, 1974). The log-normal distribution is 


n(r) — 


1 


r (2 7 r) 1 / 2 


-(lnr— lnrg) z 
2cr g 


(15) 


where r is the aerosol radius and n (r) is the number density 
of aerosols at radius r. r g and er g are the fit parameters, and 
are used to derive the effective radius, r e , and variance, u e , 
with 


r e — r g (1 + i> e ) 5 ^ 2 

(16) 

Ue = e°l - 1 . 

(17) 
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Fig. 2. Sample results from an optimization using HSRL data. In the top panels, measured reflectances, and their associated errors, are 
plotted with solid lines with respect to the view zenith angle. Model results are plotted with dashed lines. Negative view zenith angles 
indicate aircraft forward scan directions, while positive values represent aft observations. Rq is plotted on the left, while Rj is on the right. 
The bottom panels display the measurement minus model residuals. Dotted lines in these panels indicate the total error associated with that 
observation. 410mn results are in blue, 470nm in purple, 555 nm in cyan, 670nm in green, 865 nm in red and 1590 nm in magenta. 


This fitting was performed to avoid the noise sensitivity 
associated with direct computation of effective radius and 
variance (Eqs. 2.53 and 2.54 in Hansen and Travis, 1974), 
which are functions of the second, third and fourth moments 
of the size distribution, and therefore more prone to noise. 
Measurement uncertainty is the uncertainty associated with 
the fitting; observations with large uncertainty in either ef- 
fective radius or variance indicate that the aerosol size distri- 
bution is not well represented by a log-normal distribution. 
Furthermore, fitting was performed with a bimodal size dis- 
tribution (the total number concentration was the summation 
of two terms like that in Eq. 15). However, since the fine 
mode dominated the coarse mode to such a degree that the 
latter was difficult to even detect, subsequent comparisons of 
RSP and HiGEAR size distributions apply to the fine mode 
only. 

Scattering and absorption coefficients are necessary to 
compute the Angstrom exponent and single scattering 
albedo. Scattering coefficients were measured at 450, 


550 and 700 nm by a TSI model 3563 nephelometer. Data 
were corrected for ambient temperature and pressure and rel- 
ative humidities within the instrument were low, no more 
than 30%. The Particle Soot Absorption Photometers 
(PSAP) measured the absorption coefficient at 470, 530 and 
660 nm. Measurements were corrected according to the 
methodology of Virkkula et al. (2005), and represent values 
at ambient temperature and pressure. The Single Scattering 
Albedo is calculated by combining both measurements with 
the ratio ur = cr s /(cr s + a a ). This was done after interpolat- 
ing the scattering and absorbing coefficients to 532 nm. The 
nephelometer has a higher measurement frequency than the 
PSAP, so temporal averaging of the scattering coefficient was 
also required. 
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3 Results 

The primary goal of this paper was to test how observations 
from lidars such as the HSRL can be used to improve the 
analysis of RSP observations of optically thick smoke. We 
therefore performed a number of optimizations that incorpo- 
rated HSRL data, then repeated them without that data and 
compared the results. This approach was taken so that the 
variability of retrieval results across the scene can be used to 
investigate the success of the method, and so the results can 
be compared in a statistical sense. The HSRL data that were 
incorporated into the RSP analysis included both an initial 
estimate of aerosol number concentration (derived from the 
HSRL optical depth at 532 nm and the initial guess of aerosol 
size and refractive index) and the aerosol vertical distribution 
defined as layers. 

Table Al in the appendix lists the time and location of 
the individual data segments. Of note is the large difference 
(between 20° and 30°) between the aircraft heading and the 
actual ground track. This large “crab” angle is due to the 
high winds experienced by the B-200 because of the prox- 
imity to the polar jet that day. Since the RSP was therefore 
unable to scan in the plane of forward motion, the spatial ex- 
tent of a data segment, when projected at the aerosol altitude, 
was larger in the cross track than the along track direction. 
This can be seen in Fig. 1, where the white polygons indi- 
cate the spatial extent of each data segment. For this rea- 
son, geometry data from the aircraft Inertial Monitoring Unit 
(IMU) was crucial. Unfortunately, we also found that there 
was occasionally some uncertainty associated with the IMU 
heading angle. A small correction to the heading angle was 
applied in order to ensure that the linear polarization angle 
remains parallel or perpendicular to the plane of scattering. 
This is what we expect for Fresnel (single interaction) po- 
larized surface reflectance, so we use 2250 nm observations 
at scattering angles away from backscatter for this purpose, 
as those observations have the largest surface contribution to 
total reflectance. 

Tables 3 and 4 contain the mean and median values for 
retrieved and derived parameters, respectively, for retrievals 
with and without HSRL data. Table 5 contains comparable 
results from the validation data. Both mean and median val- 
ues are presented, as differences between them indicate the 
possibility that values are not normally distributed. This is 
the case for several of the parameters for retrievals with- 
out HSRL data, such as the effective variance (u e ,f), fine 
mode number concentration («f), and scattering and extinc- 
tion cross-sections (er s and cr e ). Since equivalent parameters 
measured by in situ instrumentation do not indicate a non- 
normal distribution, this is the first indication that the re- 
trievals performed without the HSRL may not be successful. 

Figure 3 shows latitude indexed comparisons between 
RSP, HSRL, AATS and HiGEAR data, essentially the results 
as the aircraft flew through or above the smoke plume. Lati- 
tude was used as the reference because of the slight temporal 


Table 3. Directly retrieved aerosol parameters. The mean and me- 
dian are calculated for the set of optimizations. Uncertainties (de- 
noted ±) are the median value returned by the optimization routine. 
Ratios are the standard deviation of the set of optimization results 
divided by the median optimization uncertainty. 


Parameter 

without HSRL 

with HSRL 

median 

mean 

median 

mean 

7 Z (mf) 

1.57 

1.55 ±0.075 

1.44 

1.45 ±0.054 

l(m f ) 

0.015 

0.016 ±0.0064 

0.005 

0.005 ±0.0036 

r e,f 

0.10 

0.11 ±0.012 

0.14 

0.14± 0.018 

v e,f 

0.38 

0.32 ± 0.049 

0.23 

0.24 ± 0.048 

nf 

82.25 

61.74± 0.182 

14.27 

17.04 ± 0.111 

n c 

0.0001 

0.0001 ±0.0020 

0.0009 

0.0008 ± 0.0020 

h 

5475 

5457 ± 0.004 




difference between the B-200 and P-3 aircraft. The aircraft 
were traveling generally to the south, which in these plots is 
from left to right. HSRL data are the source of the imagery 
in each plot. The top panel is the aerosol backscatter coef- 
ficient at 532 nm (km -1 sr -1 ), which clearly shows elevated 
layers of aerosols at 3 to 4 km. While most of the backscat- 
ter in the profile occurs at the elevated level, the presence of 
backscatter below that layer indicates that aerosols are dis- 
tributed down to the surface. The dashed white line overlay 
in this plot is the top of the aerosol layer as retrieved from 
RSP data when HSRL data are not utilized. Clearly, this 
value is inconsistent, and rarely matches HSRL observations. 
Red line overlays indicate the layers derived from HSRL data 
using Eq. (14), while the orange dashed line near the bottom 
of the panel indicates the P-3 altitude (and thus the altitude 
of AATS and HiGEAR observations). 

Our validation rests on several assumptions about the ver- 
tical distribution of aerosols that we acknowledge are some- 
what contradictory. While these contradictions may physi- 
cally exist, we believe that their radiative impacts are small 
enough so that they may be neglected. These assumptions 
are: 

1. the P-3 aircraft was flying in aerosols with the same 
properties as observed from above by the RSP on the B- 
200, so that HiGEAR observations provide a valid test 
of RSP retrieval results, 

2. the P-3 is flying under a sufficiently large portion of the 
column aerosol burden that AATS optical thickness can 
be compared to RSP retrieved results, 

3. aerosols in the direct solar beam measured by the up- 
ward looking AATS on the P-3 are the same as those 
observed by the downward viewing RSP on the B-200, 
and 

4. that the temporal difference between the P-3 and B-200 
observations are minimal. 
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ARCTAS dense smoke, June 30, 2008 
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Fig. 3. RSP, HSRL, AATS and HiGEAR data are plotted with respect to latitude (aircraft were flying to the South). The top panel shows the 
HSRL aerosol backscatter coefficient at 532 nm. The dashed white line is the retrieved aerosol layer height for RSP observations that did not 
include HSRL data. Dashed red lines indicate the HSRL determined aerosol layer heights used for RSP retrievals that used HSRL data. The 
orange dashed line indicates the altitude of the P-3 aircraft, which carried the AATS and HiGEAR instruments. The second panel is the total 
aerosol optical thickness at 532 nm for RSP retrievals without HSRL data (black). RSP retrievals with HSRL data (red), AATS (blue) and 
HSRL (magenta). That color scheme is maintained for the rest of the panels, with green indicating HiGEAR data in the third, fourth and fifth 
panels (Single Scattering Albedo, Effective Radius and Effective Variance, respectively). The sixth panel, with real refractive index, contains 
RSP data alone, indicating the validation difficulty of this parameter. Finally, the bottom panel is the information content (solid lines) and 
squared retrieval error (dashed lines) for the two RSP optimization methods. 
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The HSRL data provide some indication as to the valid- 
ity of these assumptions. For example, the HSRL extinction 
to backscatter ratio at 532 nm is nearly uniform (about 50 sr) 
within the plume, indicating that the aerosol optical proper- 
ties, at least as they are expressed in that ratio, are vertically 
consistent. Furthermore, the total optical depth from HSRL 
and AATS are also quite similar (the median values differ by 
0.033 at 532 nm), which means that the AATS is observing 
a significant portion of the total aerosol column. At the so- 
lar zenith angle for our scene, the horizontal distance at the 
top of the aerosol layer between a zenith view from the P-3 
and where the direct solar beam illuminating the P-3 enters 
the layer is on the order of two kilometers. Since the spa- 
tial resolution of a single scan at aerosol layer height is 65 m 
(see Sect. 2.1), and given that aircraft crab angles enlarge the 
spatial resolution in the cross track direction (see earlier in 
this section and the white sample polygons in Fig. 1), we can 
expect that AATS observations are physically close to RSP 
and HSRL observations. 

Another assumption we make is that spheres are an appro- 
priate geometric model for scattering by aerosols. We are 
confident that this is the case, since the HSRL observes de- 
polarization ratios of about 0.05 (at 532 nm) throughout the 
column, indicating that the aerosols can be treated optically 
as spheres (Hair et al., 2008). 

Comparisons between RSP, AATS, HSRL and HiGEAR 
data are shown in the rest of the panels in Fig. 3. In each 
panel, RSP optimizations without HSRL data are in black, 
while optimizations using HSRL data are in red, and are stag- 
gered slightly with respect to latitude to aid visualization. 
AATS data are presented in blue, HSRL data in cyan, and 
HiGEAR data in green. 

While this will be covered in more detail in Sect. 4, it is 
clear in Fig. 3 that the RSP retrievals without HSRL data 
are converging to two, significantly different, types of so- 
lutions. The first solution, which is similar to what is re- 
trieved from RSP when HSRL data are used, involves low 
absorption and narrow size distribution aerosols, with real 
refractive index values less than 1.55. This agrees with re- 
sults from the HiGEAR data for single scattering albedo and 
fine mode effective variance, and to a lesser degree effective 
radius. A second solution, retrieved only when the RSP is 
not using HSRL data, are aerosols that are very absorbing, 
very small, with large effective variances and real refractive 
indices. Since the results of this second type of solution do 
not agree with HiGEAR, we assume that the optimization 
has found a false minimum in the cost function representing 
aerosols that may not be physically realistic but are optically 
a reasonable match to the RSP observations. Presumably, 
this false minima is not found when the optimization uses 
HSRL data because those optimizations are started with a 
physically realistic fine mode number concentration. How- 
ever, despite their size and refractive index differences, the 
aerosol optical thicknesses of both states are quite similar 
and match HSRL and AATS observations of spectral optical 


ARCTAS dense smoke, June 30, 2008 



Fig. 4. Here we compare spectrally dependent aerosol optical thick- 
nesses for AATS (solid) and RSP retrievals with and without HSRL 
data (dash-dot and dashed, respectively). Values are the geometric 
mean over the time and space shown in Fig. 3. Vertical lines indicate 
one standard deviation of the retrieval values in the region shown in 
Fig. 3, and have been offset with respect to the x-axis slightly to 
ease visibility. Angstrom exponents calculated within the 400 nrn- 
1000 nm range are a =2.35for RSP retrievals without HSRL data, 
a = 2.53 for RSP retrievals with HSRL data, and a = 2. 1 1 for AATS. 


depth extremely well (see Fig. 4) . Optical depth retrievals do 
therefore appear to be robust against the uncertainties caused 
by multiple minima in the measurement cost function. 

The bottom panel of Fig. 3 contains two measures of op- 
timization success. Solid lines are the Shannon information 
content, H s , in parameter space as defined in Rodgers (2000) 
as 



(j r C“‘ J + C^ 1 ) C a 


(18) 


The information content represents the decrease in uncer- 
tainty following optimization, or equivalently the reduction 
in entropy. It is directly related to the volume of uncertainty 
for a given confidence level before and after the measure- 
ment process. Uncertainty prior to optimization is expressed 
by the a priori error covariance matrix, C a , while the un- 
certainty following optimization is the error covariance ma- 
trix, C t projected into state space using the Jacobian matrix. 
The solid lines in the bottom panel of Fig. 3 are this infor- 
mation content, which is slightly larger for retrievals using 
HSRL data during optimization compared to retrievals with- 
out HSRL data. Alternatively, optimization success can be 
expressed in observation space by the final value of the cost 
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function as defined in Eq. (2) and plotted as dashed lines in 
Fig. 3. Like the information content, the final cost function 
value is slightly worse for optimizations performed without 
HSRL data, but does not indicate a convergence failure. 

Figure 4 presents the mean spectral aerosol optical thick- 
ness from the AATS instrument compared to both types 
of RSP optimization results. All three are nearly identical 
(RSP without HSRL a = 2.35, RSP with HSRL a = 2.53, and 
AATS a = 2.11), with the exception of the of the very longest 
wavelengths, where optical thickness values are extremely 
low. This further illustrates how the spectral optical thick- 
ness is not affected by the divergent states retrieved with and 
without HSRL data. 


4 Discussion 

The primary goal of this paper is an investigation of the util- 
ity of lidar data in constraining and improving aerosol prop- 
erty retrievals with the RSP. First, we want to determine if 
the retrieved aerosol properties, when using HSRL data, are 
at least representative of boreal forest fire smoke. This ap- 
pears to be the case. Comparisons with the ground sun pho- 
tometer derived climatologies in Dubovik et al. (2002) show 
that our results are consistent with the boreal forest biomass 
burning class. Of course, this was also the source of the ini- 
tial optimization values, so this comparison is only mean- 
ingful in the sense that it shows the results did not stray far 
from their original values. Unfortunately, there are few ob- 
servations, other than AERONET, of boreal forest fire smoke 
that simultaneously retrieve refractive index, size and con- 
centration. One additional source is the three wavelength 
lidar observations of Muller et al. (2005), who studied bo- 
real forest fire smoke transported long distances from Canada 
to Germany. Complex refractive index and single scatter- 
ing albedo both are consistent with our results. However, 
Muller et al. (2005) finds significantly larger aerosols, with 
an effective radius of 0.36 ± 0.05 pm for the entire size dis- 
tribution, compared to our 0. 14 ± 0.02 pm for the fine mode 
and 0.22 pm for the entire size distribution for retrievals that 
use HSRL data as a constraint. One explanation for this dif- 
ference is the growth associated with smoke aerosol aging 
(Reid et al., 1998; Muller et al., 2007), as our aerosols are 
much younger than those transported across the Atlantic to 
Europe. While few measurements of refractive index ex- 
ist that are more accurate than the uncertainties in our re- 
trievals, there are many measurements of single scattering 
albedo, which is closely related to the imaginary refractive 
index. Our value of 0.96 (at 532 nm) is somewhat high com- 
pared to tables in the review by Reid et al. (2005a), indi- 
cating a low level of absorption for smoke aerosols. Boreal 
forest fires are known to produce less absorbing aerosols, and 
indeed the “likely optical properties” for “Temperate/Boreal 
Forest Aged” aerosols in Reid et al. (2005a) agree within un- 
certainties, with a value of 0.915 ± 0.05. Boreal forest single 


ARCTAS June 30, 2008 retrieval comparison 




Fig. 5. These plots show results from a successful optimization 
utilizing HSRL data (dash-dot) compared to an unsuccessful opti- 
mization performed without HSRL data (dashed). This segment is 
the third from the left in Fig. 3, and RSP observed reflectances are 
plotted with solid lines. For clarity, the 470, 555 and 670 nm bands 
are omitted, although they show a similar situation: optimization 
retrieved reflectances are quite similar for both cases. 


scattering albedo at 440 nm is 0.94 ± 0.02 for AERONET ob- 
servations from Dubovik et al. (2002), while Eck et al. (2009) 
found average values at 440 nm of about 0.96. 

RSP results are compared to in situ observations of the 
smoke plume in Table 5, and Figs. 3, 4, 5, and 6. Optical 
thickness retrievals, both with and without HSRL data, com- 
pare well with HSRL column and AATS observations. Fur- 
thermore, the optical thickness spectral dependence is also 
well retrieved, as can be seen in Fig. 4. There is some dif- 
ference at the longest wavelengths, but the magnitude of 
the optical depth at those wavelengths is so small that the 
estimates agree within measurement uncertainty. The only 
surprising aspect of the aerosol optical depth comparisons 
is the large RSP retrieval uncertainties (this will be dis- 
cussed more below). “Successful” RSP retrievals (we will 
use this term to describe retrievals that used HSRL data and 
the retrievals without HSRL data that found nearly identi- 
cal parameters) of single scattering albedo (nr = 0.960) are 
higher than HiGEAR observations (which have a m equal 
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Table 4. Derived aerosol parameters. The mean and median are 
calculated for the set of optimizations. Uncertainties (denoted ±) 
are the median value returned by the optimization routine. Optical 
thickness (r), scattering, absorption and extinction cross-sections 
(cr s , °a and °e). an d single scattering albedo (nr), are all expressed 
at 532 nm. 


Parameter 

without HSRL 

with HSRL 

median 

mean 

median 

mean 


0.008 

0.019 ±0.0044 

0.033 

0.032 ±0.0196 

o a 

0.0010 

0.0010 ±0.00052 

0.0011 

0.0010 ±0.00093 

cr e 

0.009 

0.020 ±0.0048 

0.034 

0.033 ±0.0200 

T 

0.666 

0.671 ±0.3038 

0.695 

0.693 ±0.3903 

a 

2.35 

2.43 ± 0.204 

2.53 

2.52 ±0.258 

TU 

0.915 

0.917 ±0.0318 

0.960 

0.960 ±0.02 10 


to 0.922). “Unsuccessful” RSP retrievals of single scatter- 
ing albedo (those performed without the aid of HSRL data 
that are lower than “successful” RSP retrievals) are lower 
(m = 0.915) and generally have larger uncertainties. Com- 
parisons to HiGEAR observations for size distribution ap- 
pear reasonable for “successful” RSP retrievals, with slightly 
larger effective radii and comparable effective variance val- 
ues. “Unsuccessful” RSP retrievals have smaller effec- 
tive radii, and much larger effective variance values than 
HiGEAR observations. Agreement for effective radius and 
variance were within uncertainties for both “successful” and 
“unsuccessful” retrievals, although it should be noted that 
HiGEAR uncertainty values are rather large. This is most 
likely related to the deviation of the actual size distribution 
from the log normal function to which we are fitting. In 
Fig. 6 we show the size distribution for a “successful” and 
“unsuccessful” RSP retrieval (Scan 170 in Table Al) com- 
pared to the closest HiGEAR observation. While this will be 
discussed in more detail below, it is notable that all three dis- 
tributions are very similar in the 0. 1-0.6 pm range. “Unsuc- 
cessful” retrievals have an excess of small particles and the 
“successful” retrievals have a deficit of small particles com- 
pared to the HiGEAR size distribution. Finally, we should 
note that we do not have in situ observations of the real re- 
fractive index for comparison. As mentioned above, “suc- 
cessful” retrievals have refractive indexes that are similar to 
published boreal forest fire smoke values, although methods 
used to determine those published values have large uncer- 
tainty. “Unsuccessful” retrievals have much larger values 
that could possibly be associated with smoke particles con- 
taining significant amounts of Black Carbon (which would be 
consistent with the larger absorption of these retrievals). Our 
hypothesis that this is not physically the case is based upon 
the larger (complex) refractive index uncertainty of these re- 
trievals and their physically unrealistic size distribution and 
single scattering albedo when compared to HiGEAR obser- 
vations. 


Particle Size Distribution 



0.01 0.10 1.00 10.00 
Radius, microns 


Fig. 6. This plot presents the retrieved and observed size distribu- 
tions for the segment in Fig. 5. RSP results for a failed optimization 
without HSRL have dashed lines, RSP results for a successful op- 
timization using HSRL data are represented as dash-dots, and the 
closest HiGEAR size distribution has a solid line. Because of the 
differences in sampling area, the HiGEAR data were normalized to 
the successful RSP size distribution. 


The use of lidar data from the HSRL did improve RSP re- 
trievals, but not in the way we anticipated. We had expected 
that an appropriate vertical distribution of aerosols would 
create retrievals that match in situ observations better and 
reduce their uncertainties. While retrievals utilizing HSRL 
data are reasonable, about half of the retrievals performed 
without HSRL data converged to nearly the same solution 
with equivalent uncertainties (what we call “successful” re- 
trievals above). However, the other half of retrievals with- 
out HSRL data converged to an alternate solution that, while 
physically possible, is less consistent with in situ observa- 
tions and has higher retrieval uncertainties (what we deem to 
be “unsuccessful” retrievals). This appears to happen not be- 
cause of an inappropriate vertical distribution (otherwise the 
“successful” retrievals performed without HSRL data would 
not be so similar to those with HSRL data), but because of 
the initial values used to initiate the Levenberg-Marquardt it- 
eration. Generally, the selection of inital parameter values 
is straightforward if there are clues to the type of aerosols 
present in the scene. In our case, the aerosols are obviously 
boreal forest fire smoke, so we used the climatology for this 
type from Dubovik et al. (2002). Dubovik et al. (2002) do 
provide mean values of aerosol optical thickness (from which 
number concentration can be derived if the extinction cross- 
section is known), but the range of reported values is so large 
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Table 5. Retrieved aerosol results for validation data. The mean 
and median are calculated for the set of optimizations. Uncertain- 
ties (denoted ±) are the median value returned by the optimization 
routine. Optical thickness (r), single scattering albedo (nr) and the 
backscatter to extinction ratio (5) are expressed at 532 nm. 


Parameter 

median 

mean 

AATS r 

0.617 

0.614 ±0.0051 

HSRL r 

0.650 

0.658 ±0.0200 

RSP without HSRL r 

0.666 

0.671 ±0.3038 

RSP with HSRL r 

0.695 

0.693 ±0.3903 

HiGEAR m 

0.922 

0.921 ±0.0181 

RSP without HSRL m 

0.915 

0.917 ±0.0318 

RSP with HSRL nr 

0.960 

0.960 ±0.0210 

HiGEAR r e f 

0.112 

0.1 12 ±0.0022 

RSP without HSRL 7 e f 

0.100 

0.113 ±0.0118 

RSP with HSRL r e , f 

0.142 

0.140 ±0.0180 

HiGEAR u e ,f 

0.252 

0.252 ±0.0011 

RSP without HSRL r> e f 

0.380 

0.325 ±0.0492 

RSP with HSRL u e , f 

0.230 

0.238 ±0.0476 

HSRL S 

49.4 

49.5 ±4.5 

RSP without HSRL S 

54.9 

57.2 ±8.6 

RSP with HSRL S 

66.0 

65.7 ± 12.1 


that the aerosol number concentration is unconstrained by 
climatology. This is because it is not an inherent quality of 
the aerosol type itself. Furthermore, aerosol optical thickness 
is log-normally, rather than normally, distributed (O’Neill 
et al., 2000; Knobelspiesse et al., 2004), which means there is 
the possibility of very large concentrations that are far from 
the arithmetic mean. The selection of the initial value for 
the number concentration is therefore quite difficult. As we 
see here, an incorrect selection can converge to a physically 
feasible local minimum in the cost function that is notice- 
ably different from the minima that agrees with in situ ob- 
servations. The HSRL data, then, were important because it 
was used to choose the initial number concentration, and the 
result was convergence to a much more consistent retrieval 
with lower uncertainties. 

We next want to address why the specific set of parame- 
ters retrieved in “unsuccessful” cases matched RSP observa- 
tions so well. As we can see from the dashed lines in the 
bottom panel of Fig. 3, there is no obvious difference be- 
tween the squared error for “successful” and “unsuccessful” 
retrievals. Both cases are capable of matching RSP obser- 
vations well, so in this case fit quality does not distinguish 
between “successful” and “unsuccessful” retrievals. To il- 
lustrate, Fig. 5 shows the match of model results to RSP 
observations for segment 170 (see Table Al). Retrieval re- 
sults without HSRL data were ’unsuccessful’ in this segment, 
but their simulated total and polarized reflectances (dashed 


lines) are not dramatically different than those for a “suc- 
cessful” retrieval using HSRL data (dash-dot lines; RSP data 
are solid lines). It appears that both “successful” and “un- 
successful” retrievals are similar in the observation space, Y. 
Figure 6 presents an interesting clue as to why this may be 
the case. In this plot, the retrieved size distributions for seg- 
ment 170 are shown, along with the closest HiGEAR ob- 
servation. There are considerable differences between the 
retrieved and observed size distributions for very small and 
coarse mode aerosols, but a striking similarity for radii be- 
tween 0. 1 and 0.6 pm. It appears that this size range is the 
optically relevant portion of the size distribution, and that 
the optimization technique faithfully matched those sizes, al- 
beit in two very different ways. Another possibility is that 
multiple scattering effects (recall that this is a scene with a 
large optical thickness and thus significant multiple scatter- 
ing) could mask differences between the two aerosol states. 
Figure 7 presents the single scattering total (left) and polar- 
ized (right) phase functions at 532 nm for the same segment 
with (solid) and without (dashed) HSRL data. The full phase 
function is plotted in the top row, while a zoom of the scatter- 
ing angles we observe is in the bottom row. Clearly, there is 
a difference in single scattering properties between the two 
states, and this appears to be masked by multiple scatter- 
ing (as we can see in Fig. 3). The lower single scattering 
albedo for the “unsuccessful” case must be compensating for 
the larger phase function magnitude when multiple scattering 
is considered, since the angular dependence over the range of 
observation angles is similar for both “successful” and “un- 
successful” phase functions. Indeed, the high absorption re- 
trieved in the “unsuccessful” cases can be seen as a symptom 
of an incorrect size distribution, as that absorption is needed 
to account for the larger amount of scattering by the smaller 
size distribution at the observed scattering angles. 

Differences between “successful” and “unsuccessful” re- 
trievals imply that there is a trade-off between the number 
of small particles (less than 0.1 pm) and complex refractive 
index that is weakly constrained by the observations. Both 
cases represent minima that are within the range of plausible 
aerosol properties, so success is highly dependent on initial 
value. This situation appears similar to results from an anal- 
ysis of simulated ground based estimates of polarized sky 
observations (Cairns et al., 1997). Simulated retrievals were 
performed for haze and bimodal aerosol size distributions. 
Haze effective radius was correctly retrieved, but there was 
a substantial overestimate of small particles. The retrieved 
real refractive index for this case was about 0. 1 larger than it 
should be, and the effective radius was quite different. These 
errors are similar to our ’unsuccessful’ case, especially since 
both the simulation and the retrieval size distributions were 
similar in the 0.1-1 .0 pm range (see Fig. 5b in Cairns et al., 
1997). Furthermore, differences between the single scatter- 
ing phase functions were masked by multiple scattering when 
modeling an optical depth of 0.5. 
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Fig. 7. These plots show the phase functions at 532 nm of the fine 
mode aerosols for the segment in Fig. 5. Successful optimization re- 
sults using HSRL data are shown with solid lines, while unsuccess- 
ful optimization results performed without HSRL data have dashed 
lines. Plots in the left column are the total (unpolarized) phase func- 
tion, while plots in the right column are the polarized phase func- 
tion. Figures in the top and bottom rows show the same quantities, 
although the bottom rows have been zoomed to a scattering angle 
range that was observed by the RSP in that segment. 


This illustrates a weakness in the ability to retrieve the 
number of aerosols whose size is large enough to activate 
as a cloud condensation nuclei (CCN) at higher supersatura- 
tions (roughly 0.03 to 0.07 pm). The number of aerosols in 
this size range is important if we are to understand aerosol 
indirect effects on clouds, since Dusek et al. (2006) found 
that cloud nucleation is more determined by the aerosol size 
distribution than its chemical composition. However, aerosol 
plumes with number concentrations and opacity as thick as 
we have observed may also have difficulty nucleating in an 
updraft because of the strong competition for water vapor be- 
tween particles. Sensitivity to CCN sized aerosols is a topic 
of ongoing interest that needs to be addressed for polarimet- 
ric remote sensing at lower aerosol optical depths, where 
concentrations of smaller particles may be important for the 
creation of cloud droplets. 

The large uncertainty in the retrieved aerosol optical thick- 
ness must also be discussed. Aerosol optical depth uncer- 
tainty for both “successful” and “unsuccessful” retrievals is 
typically about 50 %, which in an absolute sense has an av- 
erage value 0.39 at 532 nm for retrievals that use HSRL data. 
To find the source of this error we break down the compo- 
nents of the summation in Eq. (13), which are displayed in 


Table 6. Error components for fine mode optical thickness at 
532 nm for the segment 170 retrieval with HSRL data and solar 
zenith angle error assessed at 0.5°. These are the elements within 
the summations in Eq. (13). 



n ( 

77 (mf) 

X (m f ) 

'e.f 


n i 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

77 (mf) 

0.0000 

0.0143 

-0.0002 

-0.0379 

-0.0100 

X (mf) 

0.0000 

-0.0002 

0.0000 

0.0007 

-0.0001 

r e,{ 

0.0000 

-0.0379 

0.0007 

0.1074 

0.0272 

Ve.f 

0.0000 

-0.0100 

-0.0001 

0.0272 

0.0141 


Table 6. Uncertainties in effective radius are by far the largest 
source of optical thickness uncertainty, although uncertain- 
ties from effective variance, real refractive index and the cor- 
relation between effective radius and variance are also large 
contributors. It appears that the problems associated with a 
proper retrieval of the size distribution cascade down to the 
uncertainty analysis for optical depth. Even so, optical depth 
retrievals match in situ observations very well, and do not 
exhibit a sampling variability that is consistent with the large 
uncertainty in optical depth indicated by the analysis. The 
standard deviation of the set of HSRL assisted retrievals of 
optical thickness at 532 nm is 0.05, which is far smaller than 
the mean assessed uncertainty, 0.39. This indicates that the 
assessment of optical thickness uncertainty is too high. 

It is possible that our observation error covariance matrix, 
C r, is too large. A breakdown of the relative magnitude of 
the components of Eq. (3) shows that the error due to air- 
craft pitch angle, C p , is the largest contributer to the total 
error covariance matrix. This is based upon a rather arbitrary 
selection of 0.5° as the error in pitch angle. As described 
in Sect. 3, the Inertial Monitoring Unit (IMU) onboard the 
B-200 was not very accurate (and has since been replaced), 
so we needed to correct the heading and pitch angles. We 
did so by confirming that various features in observed data 
were aligned properly, including the polarization azimuth. 
We then chose the uncertainty in pitch and heading angle 
(0.5° each) to represent our best guess for our accuracy dur- 
ing that alignment. It is possible that we chose too large of 
an error for the pitch angle. Reducing the pitch angle er- 
ror to 0.1° significantly reduces the assessed error. If this is 
done for our sample segment (170), the norm of the error co- 
variance matrix is reduced by 32 % and the optical thickness 
error is reduced from 46.6 % to 35.1 %. Pitch angle error re- 
duction affects nearly all parameters, but the greatest impact 
is for optical thickness, as is displayed in Table 7. In that 
table, we also show the uncertainties computed without any 
geometry related errors (as will presumably be the case for 
orbital sensors). Uncertainties are further reduced, although 
optical thickness still has a relatively large 28.2 % error. An- 
other large contributer to Cj was C ag , which we included to 
account for the variability of the set of scans we averaged for 
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Table 7. Percent error for the segment 170 retrieval with HSRL data. The first row represents the original uncertainty assessment. The 
second is the uncertainty if the error covariance matrix is altered to represent a pitch angle uncertainty reduction from 0.5° to 0. 1° . The last 
row represents the uncertainty if there is no contribution to the error covariance matrix from heading or pitch error. 
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Table Al. Retrieval segment location, time, and solar and instrument geometry. Altitude is in meters above sea level. Optimizations without 
HSRL data for segments marked with a * were deemed “unsuccessful”. 


Starting Scan 

index 

UTC time 

Latitude 

Longitude 

Heading 

Track 

Solar Azimuth 

Solar Zenith 

Altitude 

160* 

0 

20:28:54 

58.040 

-104.606 

259.7 

240.0 

213.2 

38.3 

8681 

165 

1 

20:29:03 

58.036 

-104.618 

256.3 

237.0 

213.2 

38.3 

8681 

170* 

2 

20:29:11 

58.032 

-104.628 

255.1 

234.5 

213.2 

38.3 

8680 

175 

3 

20:29:20 

58.028 

-104.639 

251.8 

230.9 

213.3 

38.3 

8680 

180 

4 

20:29:28 

58.024 

-104.649 

250.0 

230.0 

213.3 

38.3 

8681 

185* 

5 

20:29:37 

58.019 

-104.660 

248.2 

226.7 

213.4 

38.3 

8681 

190 

6 

20:29:45 

58.014 

-104.669 

246.5 

225.0 

213.4 

38.3 

8681 

195* 

7 

20:29:53 

58.008 

-104.680 

246.6 

224.3 

213.4 

38.3 

8682 

200* 

8 

20:30:02 

58.003 

-104.689 

246.5 

222.3 

213.5 

38.3 

8683 

205* 

9 

20:30:10 

57.997 

-104.699 

244.4 

223.0 

213.5 

38.3 

8683 

210 

10 

20:30:19 

57.992 

-104.708 

249.4 

222.6 

213.6 

38.3 

8684 

215* 

11 

20:30:27 

57.986 

-104.718 

244.1 

222.3 

213.6 

38.3 

8683 

220* 

12 

20:30:35 

57.981 

-104.728 

243.9 

223.1 

213.6 

38.3 

8685 

225* 

13 

20:30:44 

57.975 

-104.738 

240.4 

223.4 

213.7 

38.3 

8687 

230 

14 

20:30:52 

57.969 

-104.747 

242.1 

223.3 

213.7 

38.3 

8687 

235 

15 

20:31:01 

57.964 

-104.757 

242.7 

222.1 

213.8 

38.3 

8686 

240* 

16 

20:31:09 

57.958 

-104.767 

241.8 

222.6 

213.8 

38.3 

8686 

245* 

17 

20:31:18 

57.952 

-104.777 

240.9 

222.5 

213.8 

38.3 

8687 


each segment. C ag is not as large as C p , but its reduction 
could also bring the optical depth uncertainty down to a level 
that is more consistent with what one would expect based on 
sampling uncertainty. In other scenarios, averaging over a 
greater number of scans would reduce C ag , but because of 
the large spatial variability we encountered in this plume we 
did not take this additional step. 

5 Conclusions 

In this paper, we have retrieved aerosol parameters of smoke 
from recent boreal forest fires. This is a particularly difficult 
type of aerosol to retrieve because it is both absorbing and 
unevenly distributed vertically, which is information that is 
difficult to retrieve simultaneously. We tested the utility of li- 
dar data when calculating optimal estimates of aerosol prop- 
erties using our airborne polarimeter. To do so, we did about 
twenty retrievals on data collected when our instrument, the 
RSP, was overflying a fresh smoke plume. We performed 
these retrievals with the RSP data alone, and then repeated 


them using information gathered by a lidar instrument called 
the HSRL. This information included observations of the al- 
titude of the elevated aerosol layer and the aerosol column 
optical depth. The latter was used to compute the initial 
value for number concentration during optimization, while 
the aerosol layer altitude was held fixed with the structure 
determined by the lidar. 

We retrieved aerosol properties that are characteristic of 
boreal forest fire smoke when our optimizations were con- 
strained by HSRL data. These aerosol properties compared 
well to a variety of in situ observations measured concur- 
rently within the smoke plume. Retrievals performed with- 
out HSRL data were successful (in terms of matching the 
HSRL constrained retrievals) only about half of the time. 
While “successful” retrievals without HSRL data were nearly 
identical to retrievals performed with the aid of HSRL data, 
“unsuccessful” retrievals found aerosols that were small (but 
with very wide size distributions and an excess of small par- 
ticles) with unusually large real and imaginary modes of the 
complex refractive index (the latter indicating significant ab- 
sorption). These “unsuccessful” retrievals occur because the 
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initial number concentration value is too far from the actual 
solution and the optimization became trapped in a local min- 
imum. HSRL estimates of aerosol optical depth (and thus 
number concentration) were therefore an important factor in 
retrieval success, and resulted in more consistent and stable 
results for the set of analyzed scenes. 

This work presents a strong argument that operational re- 
mote sensing of high loads of small particles with instru- 
ments such as the Aerosol Polarimetry Sensor will require 
either an independent estimate of number concentration, or 
the help of a table lookup to provide an initial value. Oth- 
erwise, incorrect values may be retrieved for refractive in- 
dex and size, although it is important to note that mean re- 
trieval single scattering albedos with and without HSRL data 
are within uncertainty ranges. In all cases the spectral opti- 
cal depth retrievals were consistent with AATS-14 measure- 
ments. Despite the substantially different microphysical re- 
trievals, this retrieved property is robust. Aerosol vertical 
distribution does not appear as important as it was for the 
smoke aerosols analyzed by Waquet et al. (2009), but this 
may be due to the large vertical extent and weak absorption 
of this boreal forest fire smoke. 

We have also demonstrated the importance of an accu- 
rate understanding of aircraft geometry for retrievals such 
as this. Our error estimate of 0.5° for the view zenith an- 
gle was the largest contributor to our error covariance matrix 
and is responsible for the very large estimate of uncertainty 
for aerosol optical thickness (although it contributed to the 
uncertainty of all retrieved parameters). 

The RSP and HSRL instruments remain on the B-200 air- 
craft despite the failure of the NASA Glory satellite launch. 
Future field campaigns including in situ measurement instru- 
mentation will have the capability to repeat the comparison 
methods presented here, and hopefully establish the utility 
of future scanning polarimeter observations for a variety of 
aerosol types. 
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